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Abstract 



Applying exact QCD sum rules for the baryon charge and energy momentum con- 

^^ servation we demonstrate that if the only degrees of freedom in nuclei were nucleons, the 

structure function of a nucleus would be the additive sum of the nucleon distributions 

I— I at the same Bjorken x — AQ^ /2{pA ■ q) < 0.5 up to very small Fermi motion corrections 

1-^ if l/2miqx is significantly less than the nucleus radius. Hence QCD implies that the 

'T^ proper quantity to reveal violation of the additivity due to presence of nonnuclconic 

Q^ degrees of freedom in nuclei is the ratio Ra{x,Q'^) — {2/A)F2a{x,Q'^)/F2d{x,Q'^). 

2^ Use of variable Xp — Q^ /2qQmp in the experimental studies instead of x leads to the 

[^ I deviation of Ra{xp,Q^) from one even if the nucleus would consist only of nucleons 

with small momenta. Implementation of QCD dynamics accounts in the case of the 

light nuclei for at least a half of the deviation of Ra{xp, Q^) from one for x < 0.55. In 

the case of heavy nuclei account of the QCD dynamics and of light cone momentum 

fraction carried by Fermi, Weizsacker, Williams equivalent photons are responsible for 

^vq ~ one half the deviation of Ra{x,Q^) from one at a; < 0.55. We argue that direct 

lO observation of large and predominantly nucleonic short-range correlations (SRC) in 

^y-\ nuclei impacts strongly on the understanding of the EMC effect for x > 0.6 posing a 

f^ serious challenge for most of the proposed models of the EMC effect. The data are 

CN consistent with a scenario in which the hadronic EMC effect reflects suppression of 

^^ rare quark-gluon configurations in nucleons belonging to SRC appears to be the only 

^ viable. The dynamic realization of this scenario is presented in which quantum fluc- 

*^^ tuations of the nucleon wave function with a; > 0.5 parton have a weaker interaction 

^' ^ with nearby nucleons, leading to suppression of such configurations in bound nucleons 

C^ and to the significant suppression of nucleon Fermi motion effects at x > 0.55 giving a 

right magnitude of the EMC effect. Implications of discussed effects for the analyses 

of the neutron structure function and nuclear parton distributions are presented. The 

directions for the future studies and challenging questions are outlined. 



1 Introduction 

Nearly thirty years ago the first observation of a difference between the nucleon and nucleus 
parton distributions was reported by the European Muon Collaboration (EMC) pLj. Hence 
such a difference is referred to as the EMC effect. The EMC measurement was followed by a 
series of experiments in the eighties which provided first information on the A-dependence 
of the EMC effect [2] and the A-dependence of antiquark distributions fH]. For a review of 
the data collected during the first decade of studies of the nuclear parton distributions see 
[Ij. Interest in the EMC effect was revitalized by the recent high precision measurements at 
T JNAF for a range of the lightest nuclei [5j . To extract the EMC effect from experimental 
data one should account for both the QCD dynamics where a hard process probes partons 
and the QED physics where equivalent photons carry noticeable fraction of the nucleus 
momentum in the case of medium and heavy nuclei. 

The parallel development is the BNL and TJNAF experimental studies which allowed 
one to observe directly the short-range correlations (SRC) in nuclei and explore their struc- 
ture, for a review see [6l[7j. When combined with other experimental measurement of hard 
nuclear phenomena these observations put very strong constraints on the interpretation of 
the origin of the EMC effect which are not satisfied by many of the proposed explanations 
of the effect. 

In the review we restrict our discussion to the region x > 0.2 where nuclear shadowing 
antishadowing phenomena are not important. For the recent review of small x leading 
twist physics and extensive list of references see [8] . 

Within the leading twist approximation of QCD nuclear structure functions depend on 
the Bjorken variable x = AQ'^ /2{q ■ pa) and Q^ as the direct consequence of the parton 
structure of the light-cone nucleus wave function. Description of structure functions in 
terms of the Bjorken variable x and Q^ accounts for the experimental observation that a 
nucleus consists of quarks and gluons. Account of the QCD dynamics removes kinematical 
effect introduced by the EM collaboration and followed by all other experimental groups 
who ignored parton structure of wave function of nuclear target by using nonpartonic 
variable Xp = Q'^/2qQmp to compare nucleus and deuteron structure functions, see Fig. 2. 

In the case of medium and heavy nuclei an additional effect becomes important: a 
nucleus at rest has an electric charge Z and related Coulomb field which is the zero com- 
ponent of electromagnetic field. Under the Lorentz transformation to the frame where the 
nucleus has a large momentum, the nucleus Coulomb field is transformed into the field 
of equivalent photons. This phenomenon is known as Fermi- Weizsacker- Williams (FWW) 
approximation |9j . Equivalent photons carry a noticeable and calculable fraction of the nu- 
cleus momentum if the parameter characterizing QED phenomena Zaem is not too small. 
We evaluate this fraction which is dominated by the coherent contribution into photon dis- 
tribution in a heavy nuclear target [lOj.The photon distribution in protons and neutrons 
was evaluated previously in |11[ [T^ , while in [13j the Q^ evolution of photon distribution 
within a nucleon resulting from the Q^ evolution of quark, antiquark nucleon parton distri- 



bution functions (pdfs) has been calculated R However the implications of photon parton 
component for the momentum sum rule - one of the important phenomena in the case of 
medium and heavy nucleus - was considered only in [TO]. 

The review is organized as following. In the section 2 we explain why and how an 
account of the QCD dynamics allows us to fix dynamic variables eliminating kinematical 
effects leading to breakdown of the nucleon additivity for the nuclear structure functions 
in the approximation when the nucleus is built of nucleons only. 

In section 3 we build and analyze standard nuclear model which assumes that nucleus 
in the infinite momentum frame (or light-cone) consists of nucleons only. We use the exact 
QCD baryon number and momentum sum rules to derive the formula for the nuclear pdfs 
which takes into account nucleon Fermi motion and show that corrections to the ratio of 
nuclear structure functions due to this effect are small for x < 0.55. Hence it will be 
convenient in our theoretical discussion to consider separately two kinematical regions: 
X < 0.5 where the answer is expressed through the average deformation of the structure 
functions of the bound nucleon and x > 0.6 where interplay of quantum fluctuations of 
nucleon properties with the nucleon Fermi motion effects and SRC sets in. 

In section 4 we evaluate photon structure function of a nucleus, Pa{x,Q'^), and the 
fraction of nucleus momentum carried by equivalent photons and extend the standard 
nuclear model to include this effect. 

In section 5 we demonstrate that two model independent effects discussed in sections 3 
and 4 explain about half of the EMC effect for x < 0.55 with the QCD effect discussed in 
section 3 dominating for A < 12 and two effects giving comparable contributions for heavy 
nuclei. The remaining effect represents the genuine "hadronic" EMC effect which does not 
exceed 5% for x ~ 0.5 and much smaller for smaller x, indicating that wave functions of 
bound and free nucleon are very close for most of the quark-gluon configurations even at 
the central nuclear densities. 

In section [6] we summarize the recent direct observations of the short-range correlations 
in nuclei which indicate that probability of such correlations is large (on the scale of 20%) 
and that nucleons in the SRCs are rather weakly deformed. 

In section [7] we discuss interrelation between SRC and EMC effect and demonstrate 
that the A-dependence of the hadronic component of the EMC effect is consistent with the 
A-dependence of the SRCs. 

In section |8] we summarize the most important constraints on the models of the EMC 
effect coming from the experimental studies of hard phenomena with nuclei and explain that 
wide classes of the current models are not consistent at least with one of these constraints. 

In section [9] we consider plausible interrelation between well established properties of 
hard processes off a free nucleon and dynamics of the hadronic EMC effect. We explain 
that the observed rapid decrease of antiquark distribution within a nucleon at x > 0.4 



^Authors did not explicitly considered the total momentum fraction, though they calculated the photon 
parton distributions for all x. 



indicates the suppression of meson field of a bound nucleon participating in hard processes. 
We argue that this suppression leads to a significant suppression of the contribution of 
SRC in the hard processes which select certain quark-gluon configuration in an interactive 
nucleon. Therefore Fermi motion effects and hence nuclear pdfs at x > 0.6 are suppressed 
as compared to the standard nuclear theory giving a right magnitude of the hadronic EMC 
effect at large x. 

In section 10 we discuss briefiy implications of discussed in the review effects on global 
fits of pdfs. 

In section 11 we present the conclusions and implications of the EMC effect related 
physics for different phenomena. 

2 Account of the QCD dynamics fixes the kinematic break- 
down of the nucleon additivity for nuclear pdfs 

It was proved long ago that in the limit Q^ — )■ oo, where Q^ is the square of four- momentum 
transferred by electrons to the target, but fixed Bjorken xt, the fraction of nucleus mo- 
mentum carried by interacting parton: 

XT = QV^iq-PT), (1) 

the structure functions of any target T are given by the convolution of cross section for 
hard probe scattering off individual partons with parton distribution within this target: 

The Bjorken scaling for the structure functions and Q^ evolution follow from the asymp- 
totic freedom in QCD and is implemented within the parton model accounting for the Q^ 
evolution, the light-cone quantization of QCD and the Wilson operator product expansion. 
Normalization of the light-cone wave function (WF) of a target is derived based on the 
evaluation of the matrix elements of the exactly conserved currents between the WFs of 
the target in initial and final states at zero momentum transfer (electromagnetic current, 
charmed, bottom currents). The parton model generalized to include the Q^ evolution is 
the basis for QCD physics of hard processes, for the searches of new particles, etc. 

Thus the appropriate variables to describe the process e+A — t- e'+X in terms of nuclear 
pdfs are the Bjorken variable xa (= xt) and Q^. To simplify formulae it is convenient to 
rescale Bjorken xa by the factor A: 

x/A = Q^/{2q-pA). (3) 

So x/A is the fraction of the total nucleus momentum carried by the interacting parton 
and < x/A < 1. To investigate nuclear effects which can be interpreted within QCD the 



EM Collaboration introduced the ratio: 

Ra = {2/A)F2Aix, Q^)/F2d{x, Q2). (4) 

It is important to realize that the EMC ratio as defined in Eqj4] is presented in the 
experimental papers as a function of the variable Xp = Q^ /2qQVfip [1] (also it was normal- 
ized to the cross section of electron ~ deuteron scattering rather than to the sum of the 
electron - proton and electron-neutron cross sections). Such variable is convenient for the 
experimental studies as one can compare the cross sections of the lepton - nucleus scatter- 
ing for the same kinematics of the incident and final lepton. Variable Xp differs from the 
parton model variable x and therefore within the parametrizations of structure functions 
in terms of Xp Bjorken scaling and well established universality of hard processes obtain 
unconventional form which is in variance with the physical intuition. In particular, the 
ratio Ra{xp, Q^) defined similar to Eq. mjis different from one in the kinematics where 7* 
scatters off bound nucleon carrying xa = 1/^ fraction of nucleus momentum, i.e. off a 
bound nucleon at rest in the nucleus target rest frame. It has easily calculable dependence 
on A. Thus the baryon charge and momentum sum rule are violated within a model where 
nucleus consists of nucleons only etc. Moreover nuclear pdfs presented as a function of Xp 
cannot be used directly without the additional correcting factor to predict hard phenomena 
in pA and AA collisions where as the consequence of the QCD factorization theorem the 
hard cross sections are controlled by the parton distributions over the Bjorken x fraction 
of the colliding nucleus energy carried by the interacting parton. 

We will demonstrate in the next section that the use of the Bjorken x instead of Xp 
explains at least a half of the EMC effect for the lightest nuclei {A < 12) for x < 0.5 
pUj . see Fig J2] below. (Note however that the overall EMC effect is rather small in this 
kinematical domain and there is a tendency for a smaller EMC effect in the Jlab data than 
in the SLAC data, though all data are consistent within the normalization uncertainties.) 
For medium and heavy nuclei a comparable contribution to the EMC ratio originates from 
another model independent effect - presence of the equivalent photon field of nuclei (section 
HI). Together these two effects explain approximately one half the EMC effect for x < 0.5 
for large A. 

3 Standard model for the structure functions of nuclei 

In order to discuss to what extent the EMC effect signals presence of new physics in the 
nuclear structure we need to establish expectations of the textbook nuclear model in which 
nuclei are build only of nucleons which have the same internal structure as free nucleons. 
We will refer to this approximation as the standard nuclear model. (Actually the standard 
model should be extended to also include equivalent photons as nonnucleon degrees of 
freedom. This will be done in the next section.) 

At the first step of our analysis the Fermi motion of nucleons can be neglected. Let us 
consider a nucleus moving with a large momentum P and a parton which carries a fraction 



XA of the nucleus momentum. Each nucleon carries fraction P/A of the nucleus momentum 
since the Fermi motion of the nucleons is neglected. Hence this parton carries a fraction 
Ax A of the nucleon momentum. As a result we find 

f\[xA, Q^) = Zf^iAxA, Q^) + NfiiAxA, g'), (5) 

where xa = Q^/2((/ • P)- In the case of deep inelastic scattering necessity to use xa follows 
from the fundamental property of parton model that 7* interacts with an individual parton. 
Note that the exact QCD sum rules for the baryon charge and momentum conservation 
were derived using EqJ5]or by direct calculation of matrix elements of conserved currents: 
the baryon charge and the energy- momentum tensor. Since in this approximation /^ is 
equal zero for xa > ^/A it is convenient to rescale 

XA — ^ AxA = X, (6) 

so that X is now changing between and A and consider the ratio 

^' '^ ' Zf^{x,Q^) + Nf^{x,Q^y 

If there are no nuclear effects, 

R^^{x,Q^) = l. (8) 

In the nonrelativistic approximation to the nucleon Fermi motion, the nucleus cm. 
motion is separated from the inner motion of nucleons. So the light- cone WF of a nucleus, 
tpA = expi(p^ • x)ipint, can be conveniently calculated in terms of equal time WF of the 
nonrelativistic nuclear theory - ipmt- Here pA is nucleus momentum and x the coordinate 
of the nucleus center of mass. Thus the Lorentz boost in this approximation is trivial ~ 
it is reduced to the transformation of the plane wave. This property allows us to use the 
experience of the low energy nuclear physics to evaluate some hard QCD phenomena. 

There exists a variety of the exact sum rules for valence quark distribution which 
follow from the Ward identities and existence of exactly conserved currents: baryon charge, 
isotopic charge etc. The typical valence quark sum rule for the baryon charge and the 

momentum sum rule are: 

A 

dxVAix,Q^) = B, (9) 



r-A 

/ dx[xVA{x,Q'^) + xSA{x,Q^) + xGA{x,Q'^)]=A. (10) 

Jo 

Here xVa{x,Q^) is the density of valence quark distribution, xaSa{xa,Q'^) is the density 
of nucleus sea quark distribution, xGa{x, Q"^) is the density of nucleus gluon distribution. 



B is nucleus baryon charge. An additional factor of A in Eq. 10 reflects our choice of scale 



for X - the fast momentum is measured in units of pa/^ that is average nucleon momentum 
(in the model where nucleus is build only of nucleons in the fast frame) . 

In this section we consider nucleus as a system of nucleons and use the above sum rules 
to derive approximate formulae for the nuclear pdfs accounting for relativistic corrections 
due to nucleon Fermi motion. In the impulse approximation nuclear pdfs are described by 
the convolution formulae : 



/A(x,Q') 



da 



a 



[fl{x/a,Q^)pl{a) + mx/a,Q^)p^Ac^)] . 



(11) 



Here p^"'{a) is the proton (neutron) light-cone densities of the nucleus, a/ A is the fraction 
of nucleus momentum carried by interacting nucleon, Z and N are the numbers of protons 
and neutrons. The evaluation of the matrix element of the conserved currents (electromag- 
netic, isotopic, baryon charge) - {A\ J^{t) \A) at the zero momentum transfer, t between 
nucleus states gives 

da 



A 



//( 



a) 



Z{N). 
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Evaluation of the matrix element of the energy-momentum tensor gives 



A 



[PA(a) + P\{a)]da = 1. 



(12) 



(13) 



Since the nucleus is a nonrelativistic system, the light-cone density is concentrated near 
a = 1. Hence it is legitimate to evaluate the contribution of the Fermi motion for x < 0.7 
by decomposing integrand of EqfTT] in the Taylor series in powers of a — 1 and account for 
the first three terms in the expansion. Thus the contribution of the Fermi motion effect to 
the ratio defined in Eq. [7] is [HI [15] : 



Ra{x,Q^ 



1 

A 



p^Xia) 



da xfj^'{x,Q^) 



+ 



a 



ff{x,Q^) 



1 



Pa {a)da/A 







+ - 



x/f'(x,g^) + f/f"(x,Q^)l 



/f(x,Q2) 



A 



PA{a){l-a) —. 
" 



(14) 



Here / = {f^ + fj')/2 and we assumed equal number of protons and neutrons. In the 
first term in the above equations the nucleon structure functions are cancelled in the ratio 
since the variable x is independent of the atomic number the resulting factor is equal to 
one due to the normalization condition Eq]12[ 

It is worth emphasizing that the discussed decomposition becomes inapplicable at large 
X since it diverges for x ~ 1. In this region Ra{x , Q'^) / Ra' {x , Q"^) ~ {Ta)/{Ta/) which is 
qualitatively different from Eq|14[ Here {Ta) is the average kinetic energy. 



It is convenient to introduce 



VA = l- Pa ia)da/A = {1 - a) , (15) 

Jo 

the fraction of the hght-cone momentum which is carried by nuclear constituents other than 
nucleons. It is equal to zero in the many nucleon approximation, but we keep this term 
in anticipation of the discussion of the Coulomb field effect below. In the nonrelativistic 
limit k'^/mj^, eA/rriN ^ 1 (cA is the nuclear binding energy per nucleon) the nuclear factor 
in the third term of Eq, 14 ((a — 1)^) = A~^ J^ p^{a){l — a)^^ can be calculated using 
nonrelativistic approximation for a ~ 1 + kz/rriN where k is nucleon three momentum 
within the nucleus target. (The formula a ~ 1 + h^/mN + 0{k^/m'j^) follows from the 
analysis of Feynman diagrams as the consequence of the large mass of a nucleon as compared 
to the nucleon Fermi motion momenta.) This leads to 



3^ + ^^^'/"^^'^^/"^^) = S- 



((a - 1)^) = -^ + 0{kym%, e\lm'^) = — ^. (16) 



This allows us to rewrite Eq. [T4|as 



In the case of fj^{x, Q^) oc (1 — x)", Eq. 
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can be written as 



Ra{x,Q^) = 1 + ^^^ ^ '---^-r^AUx {l-x). 18 

(1 — xy iruN 

We will present numerical results in section [5} Here we just note that it follows from EqfTsl 
that for X < 0.5 effects of Fermi motion are small and under control as the consequence of 
the account for the exact conservation laws. The Fermi motion gives zero contribution to 
Ra at X = 2/(n + 1), a small negative contribution at x < 2/(n + 1) and leads to a rapid 
growth of Ra at X > 2/(n + 1). 

At the same time the experimental data on nuclear pdfs are usually presented in terms 
of the variable 

xp = QV(2"ip'?o), (19) 

which is different from the fraction of target momentum carried by interacting parton and 
therefore requires rerwriting formulae for parton model and Q^ evolution. This variable 
depends on the atomic number of a target even if the nucleon Fermi motion effects are 
neglected. 

Xp/x = Anip/mA = (1 + (e^ - (m„ - mp)N/A)/mp) = 1 + r^,. (20) 



Thus the structure functions of interacting nucleons in the nominator and denominator de- 
pend on different arguments. Thus the use of variable Xp instead of x (and the expectation 
of absence of the nuclear effects in such an approximation) violates important condition of 
standard nuclear model that nuclear pdfs are the sum of nucleon pdfs if the Fermi motion 
effects are neglected and introduces artificial dependence of Ra on atomic number. In 
particular RAix) = 1 corresponds to 

Rixp) = fi{x{l + r.,))/fj,{x) « 1 - r^n-^, (21) 

1 — x 

where at the last step we took fj^{x) oc (1 — 3;)"rl 

The nucleon Fermi motion correction can be easily included as well. We will demon- 
strate in section [5] that Eq|21| explains at least half of the EMC effect for x < 0.5 for light 
nuclei A < 12 - i.e. the EMC effect in this case to a large extent is due to dependence of 
Xp on the atomic number. With a further increase of A the deviation of Ra{x) from one 
grows. We will explain in the next section that a large fraction of this growth is due to the 
presence of the Coulomb field in nuclei which strongly grows with A. 

Comment. A popular expectation in the low energy nuclear physics is that one can 
account for the effects of relativistic nucleon Fermi motion assuming that the vertex func- 
tions in the Feynman diagrams with a virtual nucleon coincide with Schrodinger WFs of a 
nucleus. This model has been applied in several papers to explain the EMC effect formu- 
lated in terms of non-parton model variable Xp and without subtraction of the contribution 
of equivalent photons. In the first papers [16l[T7] the baryon sum rule, i.e., the probability 
conservation was violated. Later the baryon charge sum rule has been taken into account 
following the prescription of [18j. This non-partonic description does not allow one to 
account consistently for the probability conservation i.e. to satisfy simultaneously both 
baryon and momentum sum rules for the hadronic part of the EMC effect. This is because 
space-time evolution of hard processes in QCD is described by the light-front dynamics. 

The baryon sum rule in this model is derived by calculating matrix element of baryon 
charge exactly i.e. by modifying nonrelativistic normalization of the nucleus WF. n (Note 
that normalization of wave function natural for the nonrelativistic physics implies correction 
to the total cross section of DIS off deuteron of the same magnitude as the Glauber shadow- 
ing correction which violates probability conservation ~ the so called West correction|19j'l. 
In this approach there is no symmetry between distributions over fractions of nucleus mo- 
mentum carried by interacting nucleon and nucleon spectators. As a result in this model 
significant light-cone fraction of the nucleus momentum is carried by the non-nucleonic 
degrees of freedom [15] : 

r]A = eA/rriN + Ta/^itin, (22) 



■^Note that to simplify the expressions we took here m„ 



the final expressions the A/D ratio we take into account the difference between Xp and xd- 

■^For the detailed discussion of the early theoretical studies of these issues and references see |15) . 



For the realistic nuclear WFs and A ~ 200, rjA ~ 4% which leads to a large suppression of 
-Ryi(x ~ 0.5) which is in addition to the standard model effect of the QCD dynamics choice 
of the Bjorken variable and the effect of due to presence of the equivalent photons discussed 
in the next section. After the two model independent effects we discuss in the review 
above are included, the virtual nucleon model leads to significantly larger suppression of 
the nuclear structure functions at x < 0.5 than observed experimentally. 

4 Photon distribution in nucleons and nuclei. 

4.1 General framework 

The long range Coulomb field of nucleus at rest is the zero component of the electro- 
magnetic field. Lorentz transformation of nucleus Coulomb field from rest frame to the 
frame where nucleus is rapid unambiguously leads to the field of Fermi, Weizsacker and 
Williams equivalent photons [9]. Thus the light-cone WF of nucleus contains photons as 
constituents carrying a fraction of the nucleus momentum. Collision of equivalent photons 
of one nucleus with nucleon or nucleus beams produces variety of hard and soft high energy 
processes in the ultraperipheral processes at the LHC [20] . Ultraperipheral processes were 
observed at RHIC, see review in j20j . 

For this paper a proper example of the ultraperipheral processes is the diffractive pro- 
duction of massive lepton pairs 7* + A — t- L^ + L~ + A. Large Q^ and/or large mass of 
the lepton pair L^L^ guarantees dominance of the leading twist term and allows one to 
define photon distribution in a nucleus, PAix^Q"^). 

A nucleus is characterized by quark, gluon, photon distributions within a nucleus. To 
suppress particle production by hard probe from vacuum the gauge condition A'^ = is 
chosen, where Afj_ is the operator of photon field. In this gauge photon distribution has the 
form familiar from QED, cf. Fig. [T] The photon parton distribution can be written as the 
matrix element of the product of the operators cf . [21j : 

1 f"" 

Pa{x, Q2) = -(27rxp+)-i / dy- exp {-i{x/A)p+y-) 
71" J-00 

■l/2Y,{A I [f;(0, y- , 0), F'^'+IO)] | A)a+=,. (23) 

Here i^^,+ is the operator of the strengths of the photon field with transverse component 
fi, and x/A = Q"^ /2{pAq) is the fraction of nucleus momentum carried by a parton. 

To simplify calculations it is convenient to represent Pa as the sum of two contributions 
Pa = P^"*^' + pcoherent ^ First term accounts for the target excitations by virtual photon - 
we will refer to it as the inelastic term. The second term is the contribution of equivalent 
photons. We will refer to it as the coherent term. 

Small value of Uem guarantees that in the kinematic domain achievable at accelerators 
aemlii((5^/-^n(7_D)/^^ ^ ^ ^"^^ therefore all effects of this order can be neglected. In 
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Figure 1: Diagram for the interaction of photon of the nucleus with a hard probe h. 

particular, in this kinematic domain the effects of running of aem{Q'^) due to combined 
QED and QCD effects are negligible. Moreover it is legitimate to neglect in this kinematics 
by the Q^ evolution of the amplitude of the scattering of a hard probe h off the photon which 
leads to tiny corrections. There are two practical consequences. The Q^ evolution of photon 
distribution in a target T - Px{x,Q'^) goes in one direction: quarks, antiquarks radiate 
photons in the process of evolution but a photon does not radiate them. This evolution is 
accounted for in the evolution equation through the Q^ dependence of valence, sea quark 
and gluon densities. Second consequence of the smallness of the electromagnetic constant 
is that coherent contribution is not evolving with Q^ at Q^ achievable experimentally. 
Above we defined p^"^' exactly in the same way as other parton densities, cf. Eq. 



23 



This definition allows to calculate P^"^' in terms of quark, antiquark nuclear pdfs by eval- 
uating corresponding Feynman diagrams. 



xPT\x, Q')/A = ^^ / de, / — ,,2^n2 'I ^2 -^2aK ^')M- (24) 

Within the leading Oe.m. log Q^ approximation upper limit of integration over photon trans- 
verse momentum k]:^^^ < Q^. Here v = 2{pq)/A = Q'^/x, kt is the transverse momentum 
of the photon, and F2a{i^, Q"^)/^ is the virtual Compton amplitude normalized per nucleon. 
In the above formulae we neglected the small contribution of the longitudinally polarized 
photons (F^). 

The presence of equivalent photons in the nuclear WF leads to the violation of the 



intuitive prediction we discussed in sections 2 3 - R^{x, Q ) = 1, for x < 0.5 



^, . ^2^ ^ Zf^{x/il-v,),Q^) + Nfi{x/{l-r,,),Q') 
^ Zfl{x,Q^) + Nfl{x,Q^) 

Here r]^f{A) is the fraction of nucleus momentum carried by equivalent photons calculated 
in [lOj. Account of the presence of tj^^^A) leads to a definite dependence of i?'^ on Z and 
A. For example this effect is as large as the effect of accounting for the difference between 
x and Xp for heavy nuclei while it is negligible for Z=2. 

11 



4.2 The coherent and incoherent contributions 

First we will calculate coherent contribution to parton nucleus distribution which domi- 
nates the photon distribution in a nucleus. This contribution to Pa{x, Q^) arises from the 
interaction of a hard probe with a photon coherently emitted by the target so that the tar- 
get remains intact, cf. FigjTj The coherent contribution to the photon structure function 
is unambiguously calculable in terms of the electromagnetic form factors of the nucleus 
target. For a proton target coherent term has been calculated in [11]. Inelastic contribu- 
tion for a nucleon target (where nucleon is excited in the final state) was calculated in the 
perturbative QCD model where quarks and gluons are generated via evolution starting at 
very small Q^ \12}- Comparison of two contributions was performed at ^13j where it was 
found that the elastic contribution is much more important in the proton case in a wide 
range of the x, Q^. In the neutron case incoherent contribution is more important. 

An important aspect of the photon contribution - its implications for the momentum 
sum rule was not discussed in these papers and effects of the presence of photons in addition 
to quarks and gluons at the normalization scale are still not included in the pdf studies by 
the groups which analyze the hard processes. We calculate the field of equivalent photons 
through the evaluation of Feynman diagram in Fig{T} Calculations are simplified in our 
case since the nucleus is heavy so the static approximation should be sufficiently accurate. 
In the static approximation zero component of photon momentum in the nucleus rest frame 
is negligible: /cq = k'^ /2mA ~ 0. So 

X = A{ko - k^)/MA « -k^/niN. (26) 

Second simplification arises from the observation that four- vector k can be decomposed 
over directions defined by external momenta: A;^ = ap^ + bq^^ + kf. Here p is four momentum 
of target nucleus and q is four momentum of virtual photon (external hard probe) and 
(ph) = (gki) = 0. In the essential region: p^Afj_^x ~ {kt^/a)A^^\, cf. [22]. Account of this 
property leads to the generalization of the Fermi - Weizsacker -Williams expression for the 
spectrum of the equivalent photons: 

jcoherentr^ ^2\ _ (^em Z f jr2i2^A\^t +^ f^N) 



t: A J [k^+x'^mj^Y 

Here FA{t) is the observed electric form factor of the nucleus which is equal to the product 
of the nucleus body form factor and the proton electric form factor Fp{t). Hence the 
characteristic values of k"^ are small. 

In our estimates we choose Fa in the exponential form, i. e. 

Fa(A;2 + x^m%) = exp{-r\{kf + x'^m%)/6). (28) 

Here r^ is the experimentally measured RMS nuclear radius and mjv is the nucleon mass. 
Such a form allows one to perform the analytic integration over the transverse momenta 
of the photons. 
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It is useful to compare different electromagnetic contributions. The most important 
one is the contribution of equivalent photons in which the fields of individual protons add 
coherently. This coherence leads to a larger momentum fraction carried by the photon 



field in nuclei as compared to that carried by individual free protons. Eq. [27] can be used 
to evaluate the coherent contribution of the photons to the momentum sum rule (here we 
already subtracted contribution of individual protons to the photon parton density - see 
discussion below, since this quantity enters into description of the nuclear effects we discuss 
in the paper): 

/■I 2 Z"^ 1 

r\^= \ dxxPAix,Q'^) = aem—7=—. ■ (29) 

Jo V3vr A mNVA 

To separate nuclear effects one need to subtract the contribution of Coulomb fields of 
protons. This is achieved by taking into account contribution of the incoherent break up of 
the nucleus as well as production of hadrons (inelastic processes) . Sum of the two effects can 
be calculated in the closure approximation where cross section is described by the sum of 
two diagrams presented in [23] and results in replacing Z'^F\{t) — > Z F'^{t)+ Z {Z —\)F\{t) 



in Eq. 27 leading to 



/■I , ., 2 Z(Z-l) \ Z 

rj^ = dxxPA{x,Q ) = aem—?= -. ^ T^p (^0) 

Jo V3vr A mNVA A 



The second term in Eq 30 is due to the Coulomb field of the individual protons and it 
is included in the structure functions of the proton. Hence it should not be included in the 
calculation of the nuclear effects. The contribution of magnetic form factors of neutron 
and proton is negligible because it is concentrated at larger momentum transfers than the 
scale of the nuclear phenomena. 



Taking va from the compilation of ref . |24j we find for the nuclear term in Eq 30 

r]^{^He) = .03%; r]^{^'^C) = .11%; //^(^^A/) = .21%; 7/^(^6 ^e) = .35%; r]^{^^'^ Au) = .65%. 

(31) 
To evaluate the impact of the presence of the photon component for the nuclear structure 



functions we can use Eq 18 and take tja = ri^{A). The small x nuclear shadowing effects 
modify this approximation, however they are negligible for x > 0.2 range we are interested 
in. Since deviations from the additivity are small the effect of the presence of the photons 
and other "hadronic" effects can be treated as contributing additively to the deviation 
of the EMC ratio from one. Larger numbers given in ^lOj are because of misprint in the 
formulae for coherent contribution and due to inaccurate numerical calculation. 

4.3 Equivalent photons and modification of the initial condition for the 
Q^ evolution of nucleus and nucleon pdfs 

Major impact of the presence of the photon constituents in a nucleus is through the change 
of the form of the energy-momentum conservation and violation of isotopic invariance in 
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the incoherent contribution of protons and neutrons. Total fraction of nucleus momentum 
carried by QCD partons instead of 1 becomes: 1 — i]^{A) (in this fraction the photon fields 
due to individual nucleons are included). 

The presence of the photon component in the nuclear light-cone WF leads to the certain 
modification of Q^ evolution of parton densities. Most important is the change of the 
momentum sum rule due to necessity to account for the fraction of nucleus momentum 
carried by photons: 

[{1/A){xVa{x, g2) + xSa{x, Q^) + xGa{x, Q^))] dx = 

1 - / [{1/A)xPa{x, Q2)] dx = l- r]^{A). (32) 

Jo 

The light cone momentum carried by equivalent photons is compensated by the loss of 
the momentum by the nucleons. Since the nucleus Coulomb field is the collective field 
dominated by the contribution of large impact parameters, the reduction of the light cone 
fraction is experienced by nucleus as a whole (i.e. by both protons and neutrons.) Hence 
the shift is approximately equal for protons and neutrons. 

Although r]^/{A) is small its role is enhanced by the rapid decrease of nuclear pdfs with 
x increase. 

In the case of a nucleon target one should also account for the change of the form of 
the momentum conservation and much larger contribution of photons for the proton target 
than for the neutron target leading to an isospin violating effect for pdfs of protons and 
neutrons and hence for the nuclear targets. 

5 Implications for the EMC effect for nuclear pdfs at a: < 0.5 

Here we study to what extent two effects we discussed above - the proper definition of x 
which preserves the momentum sum rule, and the account of the collective Coulomb field 
of the nucleus (of the momentum carried by equivalent photons) - which we refer below as 
the standard model of the nucleus - explain the data at x < 0.5 were Fermi motion is a 
small correction. 

In the standard model we can use Eqs. 18|21 to obtain 



4(x,) = f{{xil + r^+rj,))/fUx) = l-(^^+^7)^Y^ + "''^"i"_^J)l '^ -^^ (33) 
Since the EMC ratio is experimentally defined relative to the deuteron one needs to sub- 



stitute rx{A) by r^iD). The results of calculations using Eq 33 with Ta from [25] are 



presented in Fig. 2. One can see that Eq 33 without including Fermi motion effects ex- 



plain about a half of the EMC effect for x < 0.65. The Fermi motion changes the results 
for X > 0.55 indicating that genuine nuclear effects beyond the standard nuclear model 
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Figure 2: Contribution to the EMC ratio of the nucleus field of equivalent photons and 



effect of proper definition of x calculated as described by the second term in Eq 33 (solid 
curves) . The dashed curves include also the effect of the Fermi motion estimated using Eq. 



33 The data are from [5] - blue circles and [26] - red squires. Open circles correspond to 



W <2 GeV. 



strongly increase with increasing x. For the case of carbon the dominant contribution 
Ra{xp) — 1 comes from the r^ term, while for A ~ 200 the rx and the equivalent photon 



contributions are about equal Q Eq. 33 leads to a very small EMC effect for the ^He case 
namely for the ratio -F23_ffe(^p)/[-^22_ff(3;p) + F2p{xpy\ due to a very small energy binding 
per nucleon in ^He, which is consistent with the lack of the x dependence of this ratio as 
reported in [27]. Since the overall effect is very small it is difficult to compare the residual 
hadronic effect (that is the effect left after taking into account the equivalent photon effect 
and the use of the proper scaling variable) with the theoretical ideas due to the normal- 
ization uncertainties. However an effect on the level of 1/4 of the effect for heavy nuclei 
(the A-dependence expected from the contribution of the short-range correlations - see 
discussion below) - cannot be excluded.). 

33 also leads to a very similar contribution to the EMC ratio for carbon and He 



Eq. 



since the energy bindings per nucleon are very close and the difference of r/^ is only 0. 
which is consistent with the data within the normalization uncertainties. 

Overall inspection of Figj2] indicates that the contribution to EMC effect due to the 
modification of the quark distribution in nucleons is significant only for x > 0.5 - which 
we will refer to as hadronic EMC effect. 



''The data for heavy nuclei [26] were corrected for the difference of the number of protons and neutrons. 
However exact value of this correction requires further analysis since the EMC effect for the deuteron was 
probably overestimated in these analyses - see discussion below. 
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6 Recent progress in the studies of SRCs in nuclei 

Before discussing the A-dependence of the the hadronic EMC and its dynamical origin 
we need to discuss the information about the short-range correlations in nuclei. A very 
significant progress in this field has been reached in the last few years. As a result it has 
become possible to perform nearly model dependent tests of the role of the SRC in the 
EMC effect and also strongly constrain the models of the EMC effect. 

Hard interactions resolve local structure of the nucleus. Hence they resolve collective 
low energy scale degrees of freedom (effective interactions) and probe local structure of 
the nucleus. Since the nuclear forces are short range on the nucleus scale (r^N < "^ffn) 
the properties of a nucleon in the nuclear media are determined by its local surrounding. 
Intuitively it is clear that closer the nucleons get together, the stronger their polariza- 
tion/deformation is. Hence the configurations where nucleons nearly overlap appear to be 
a natural candidate for the large hadronic EMC effect. Therefore we need to review briefly 
the recent progress in the studies of SRCs. 

Singular nature of the NN interaction in coordinate space at small inter nucleon dis- 
tances/large momenta leads to the universal structure of SRCs and to the prediction of the 
scaling of the ratios of the cross sections of a; > 1 scattering at sufficiently large Q^ > 2GeV^ 
|14] . In particular for 1 + kp/raN < x < 2: 

Here 02 (^) has the meaning of the relative probability of the two nucleon SRCs per nucleon 
in a nucleus and in the deuteron, that is, it is the ratio of spectral functions of nucleus 
to that of deuteron. Actually it includes the contribution of triple nucleon correlation, pn 
and pp correlations. 

The first evidence for such scaling of the ratios was reported in [15]. The extensive 
studies were performed in |28j using various data taken at SLAC at somewhat different 
settings which confirmed the scaling of the ratios and for the first time confirmed the 
prediction [15] of the " super" scaling of the ratios at different Q^ - the precocious scaling 
of the ratios plotted as a function of at.n. ~ the minimal a for the scattering off two nucleon 
SRC at rest (the Fermi motion of the pair practically cancels out in such a ratio) [T^] : 



«*■-=' 2^ y^ W • ' ^^'^ 

where W'^ = 4m^ + 4goW'Af — Q^- 

The experiments performed at JLab allowed to explore the scaling of ratios in the same 
experiment. In |29tl30j the scaling relative to ^He was established. Very recently the results 
of the extensive study of the nucleus/deuteron ratios were reported in [3lj allowing a high 
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precision determination of the relative probability of the two nucleon SRCs in nuclei and 
the deuteron. The results of |31j are in a good agreement with the early analysis of [28], 
see Fig. 3. 

Several theoretical observations are important for interpretation of the scaling ratios: 
(a) The invariant energy of the produced system for the interaction off the deuteron is 
large on the scale of nuclear phenomena but small as compared to the scale characteristic 
of hadronic phenomena: W — m2jj < 250 MeV. This kinematics obviously leads to a strong 
suppression of the production of the inelastic final states. Correspondingly, scattering off 
exotic configurations which decay into excited baryon states, A's, etc is strongly suppressed 
in the discussed kinematics, (b) The closure is valid for the final state interaction of the 
nucleons of the SRC and the residual nucleus system. Only the f.s.i. between the nucleons 
of the SRC contributes to the total (e,e') cross section j281 [6]. Since this interaction is 
the same for light and heavy nuclei it does not modify the scaling of the ratios, (c) In 
the limit of large Q^ and large energy qq but the fixed ratio x = AQ'^ /2mAqo, the cross 
section is expressed through the light - cone projection of the ground state nuclear density 
matrix, p^{a), that is the integral over all components of the interacting nucleon four 
momentum except a = p-/{mA/A.), where p- = Po — {p ■ q)/\q\- The ratio of the cross 
sections reaches a plateau at x(Q^) corresponding to the scattering off a nucleon with 
minimal momentum ~ kp indicating that the dominance of two nucleon SRC sets in just 
above the Fermi surface, (d) Observation of the precocious at.n. scaling indicates that Ra 
is equal to the ratio of the light-cone density matrices of the nucleus and deuteron. It 
also strongly indicates that SRCs of the baryon charge two are predominantly built of two 
nucleons rather than some exotic states. It is worth noting here that the correspondence 
between the light cone and nonrelativistic wave functions is pretty straightforward. Hence 
the results for the ratio of wave functions in the region where pair correlations dominate 
should be close. 

To probe directly the structure of the SRCs it is advantageous to study a decay of SRC 
after one nucleon of the SRC is removed which is described by the nuclear decay function 
|14^ 115] . In the two- nucleon SRC approximation the decay function is simply expressed 
through the density matrix as the removal of one of the nucleons of the correlation results 
in the release of the second nucleon with probability of one. A series of the experiments 
was performed at BNL and JLab which studied (p,2p), (e,e'p) reactions in the kinematics 
where a fast proton of the nucleus is knocked out (see review and references in |32^ [6]). 
In spite of very different kinematics - removal of forward moving nucleon in the ^^C(p,2p) 
case and backward moving proton in the ^^C(e,e'p) case, different probes and different 
momentum transfer — i ~ 5 GeV and Q^ = 2 GeV - the same pattern of the neutron 
emission was observed: the neutron is emitted with a probability ~ 90% in the direction 
approximately opposite to the initial proton direction with the correlation setting in very 
close to kF{C) ~ 220 MeV/c. The JLab experiment observed in the same kinematics 
both proton and neutron emission in coincidence with e'p and found the probability of the 
proton emission to be about 1/9 of the neutron probability. Hence the data confirm our 
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Figure 3: Comparison of the first determination of 02 (^) based on the analysis of the SLAC 
data EHIwith the most recent Jlab measurements 1311 • 



theoretical expectation that removal of a fast nucleon is practically always associated with 
the emission of the nucleon in the opposite direction with the SRC contribution providing 
the dominant component of the nuclear WF starting close to the Fermi momentum. The 
large pn/pp ratio also confirms the standard expectation of the nuclear physics that short- 
range interactions are much stronger in the isospin zero channel than in the isospin one 
channel and hence are much more sensitive to the pion-like exchange. Saturation of the 
probability provides an independent confirmation of the conclusion that at least up to 
momenta ~ 500 -^ 600 MeV/c, SRC predominantly consist of two nucleons. 

A word of caution is necessary here. In the standard nuclear physics approach the 
intermediate A isobars play an important role in the 1 = 1 interaction channel. Usually 
they are absorbed into the definition of the low energy NN potential. However in the case 
of hard exclusive process the A^A SRCs are resolved and may show up and be comparable 
to pp SRCs. 

The structure of SRC depends on x,Q'^ and M^. Here M^ is the mass range allowed 
for the decay of SRC. For M close to 2mN possible exotic states enter into parameters of 
effective low energy Hamiltonian. With increase of M exotic states may reveal themselves 
in the decay of SRC. Thus M—2mj\f plays the role of the resolution and increasing resolution 
allows one to investigate some properties of superdense matter which is present in the inner 
core of neutron stars. Note here that the allowed phase volume was much larger in the 
BNL A(p,2pn) experiment. So the consistency of two experiments suggests that increase of 
the resolution does not lead to a drastic change in the structure of the short range nuclear 
structure. 
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7 Short-range correlations in nuclei and the A-dependence 
of the hadronic EMC effect 

We will show below that the strength of the deformation of the bound nucleon WF as 
compared to the free nucleon WF is proportional to its off-shellness. At x < 0.7 where the 
Fermi motion effects are small it is possible to average over nucleon momenta and conclude 
that the deformation of the structure function of a bound nucleon should be proportional 
to the nucleon average kinetic energy Ta- Since Ta is dominated by the contribution of the 
short range correlations [25j in nuclei we can roughly estimate the A-dependence of Ta from 
the measurements the ratio of the high momentum components in the nuclei and in the 
deuteron - a2 {A) . This ratio is equal to the ratio of the nucleus and the deuteron inclusive 
(e,e') cross sections in the x region of 1.4 < x < 1.8 where it exhibits a Q^ independent 



plateau - (see discussion in Section ro^ and in particular Eq 34). 

This information can be used to predict the A-dependence of the hadronic EMC effect 
for X = 0.5 where nucleon Fermi motion is absent for F2N oc (1 — x)",n = 3. One can 
see from Fig. [4] that the A-dependence of the "hadronic" EMC effect for this x is indeed 
consistent with the A-dependence of 02 (^) — 1. It is worth noting here that in |33| 115] it 
was assumed that all x > 0.4 EMC effect is due to the contribution of the SRCs. (Note 
here that nucleons belonging to mean field give up to 20% contribution to Ta, leading to 
a deviation of the A-dependence of Ta and a2{A).) 

To analyze the A-dependence of the hadronic EMC effect at x > 0.5 we need to take 
into account the nucleon Fermi motion effect which rapidly increases with x for x > 0.55 
and which is oc Ta in the discussed x-range (Eq]17|) and results including Fermi motion 
- dashed curves in Fig. 2a, 2b). One can see that the difference between dashed curves 
and the data is much larger than the difference between the standard model without Fermi 
motion and the data already for x=0.6. This implies right away that the hadronic EMC 
effect is roughly oc Ta- To check it with a better accuracy we can focus on the A-dependence 
of the difference between the standard model without Fermi motion and the data, which 
we denote as A(A). In Fig. [s] it is plotted as a function of 02 (^) — 1- The data are 
consistent with A(A) oc a2{A) — 1 though more accurate data would be highly desirable. 
Thus we conclude that in the range: 0.5 < x < 0.7 the hadronic EMC effect is proportional 
to the average nucleon kinetic energy which in turn is approximately proportional to the 
probability of two nucleon correlations in nuclei. For larger x this pattern should break 
down since the expansion in powers of (1-x) breaks down. In fact for x ~ 1 we expect 
F2A/F22H oc a2{A) rather than F2a/F22}j — 1 oc a2{A). 

Theoretical and phenomenological arguments in favor of the similarity of the A-dependence 
of the EMC effect and the two nucleon SRCs were first presented in \1'6\ ITB] for x > 0.5. 
New data on the A-dependence of SRCs allowed to demonstrate proportionality of the EMC 
effect to the probability of the two nucleon correlations in nuclei without employing model 
for the A-dependence of SRCs |34j . However the authors assumed that the SRCs dominate 
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Figure 4: The solid line is the result of calculation taking into account the equivalent 
photon effect and the effect of proper definition of x. The dashed line is the contribution 
of the hadronic EMC effect due to SRCs normalized for large A with the A-dependence 
oc a2{A) — 1 with a2{A) from taken from the JLab measurement |31j . 
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Figure 5: Deviation of the EMC ratio from the standard model without nucleon Fermi 
motion for x=0.6 and x=0.7 as a function of a2{A) — 1. The normalization errors of the 
data points (not shown) are typically ~2% for ^He and Au and ~ 1% for the rest of the 
nuclei. 
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Figure 6: Estimate of the ratio of the bound and free nucleon structure functions in medium 
and heavy nuclei as a function of x. 

the EMC effect for the R^ixp) ratio for all x. As we have seen above accounting for the 
standard model effects (which in most of the range studied in \^4\ -x < 0.6 - constitutes 
a half of the difference of Ra{xp) from unity) does not change this conclusion. It is worth 
emphasizing that presence of the equivalent photons helps to ease a long standing puzzle 
- why the EMC effect increases with A for A > 12 and x > 0.5 though the probability of 
the SRCs nearly saturates. 

Thus we conclude that the EMC effect at x < 0.7 is due to a combination of several phys- 
ical phenomena. For x < 0.5 at least half of the effect is explained by the standard model 
(use of the proper x and accounting for the contribution of the equivalent photons) . Bound 
nucleon deformation effect is proportional to T^, and the nucleon Fermi motion effect nearly 
compensates the deformation effect for x < 0.65. The bound nucleon deformation effect 
for medium and heavy nuclei that is deviation of Rb{x,Q^) = F2l^'^'^ {x ^ Q"^) / F^l^^ {x , Q"^) 
from unity is ~ 4% for x ~ 0.5 and starts to grow rapidly with further increase of x: ~ 
15% for x=0.6, ~ 25% for x=0.7 (The value of 1 — i?f,(0.7) is somewhat overestimated since 
for X > 0.7 factorization of the correction to the EMC ratio as a product of the nucleon 
modification effect and the Fermi motion and the standard model effects breaks down, 
leading to overestimate of 1 — Rb{x,Q'^) for such x. A separate analysis is required which 
will be presented elsewhere.) The estimate of 1 — Rb{x,Q^) is presented in Fig. p^ based 
on the analysis of the SLAC data|26j. 

Since the bound nucleon deformation effect is small up to x ~ 0.6 and the probability 
of having a quark in a nucleon with x > 0.5 is ~ 2 • 10~^ we conclude that the EMC 
effect probes very rare deformations of bound nucleons. The deviations from the discussed 
approximation on the level of 1 -^ 2% at x ~ 0.2 -^ 0.4 cannot be excluded experimentally 
so one can only conclude that the accuracy of the approximation where nucleus consists of 
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nucleons whose structure function coincides with that for free ones (+ equivalent photons) 
is on the level of few % for x < 0.5 where the contribution of the mean field approximation 
to the nucleus WF dominates. 

Our analysis also allows to put an upper limit 

r/^(A~200) <2.0%. (36) 

on the fraction of energy carried by pion constituents of nuclei assuming that all deviation 
of Ra{x ~ 0.5, Q"^) from the standard model is due to the pion field. 

8 Constraints on the models of the EMC effect 

In this section we briefly summarize the constraints on the possible mechanism of the 
hadronic EMC effect and confront different classes of models of the EMC effect with these 
constraints. Majority of models ignore the effects described by the standard model which 
describes half of the EMC effect and hence should be modifled to describe data. So we 
restrict our discussion by consideration of constraints on the basic features of these models. 
The most important of these constraints are obtained from the following: 

• QCD dynamics of hard processes and parton structure of nucleons and nuclei leads 
to unique definition of the Bjorken x. 

• Probability conservation: the baryon and momentum sum rules. 

• Fraction of nucleus momentum is carried by equivalent photons and accurately cal- 
culated in QED. 

• The hadronic EMC effect is small for x < 0.5 and for larger x rapidly increases with 
increasing x, see Fig. 6. 

• The A-dependence of the hadronic EMC effect is rather close to the A-dependence 
of the the short-range internucleon correlations, cf. discussion in section 7. 

• High probability of short-range correlations ~ 20 -^ 25% and dominance of the nucle- 
onic degrees of freedom (> 80%) in the correlations, cf. brief review of existing data 
in section 6. 

• No enhancement of the antiquark distribution in nuclei is present for x ~ 0.1 with 
accuracy ss 1% [3]. 

• Q^-dependence of the magnetic form factor of bound nucleons with small momenta 
is very close to that of free nucleons [35] . 
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A natural treatment of the EMC effect within conventional nuclear physics approaches 
is to explicitly include the mesonic degrees of freedom in the nuclear wave function. To 
fit the data authors assumed that mesons carry r;^ ~ 4% fraction of the nucleus light 
cone momentum. A significant enhancement of the pion field at large pion momenta was 
found for example in |36]. However relating this enhancement to the value of rj^^ is not 
straightforward as it involves a Lorentz boost and switching from the equal time to light 
cone quantization. The models which describe nucleus wave function by a vertex function 
with interacting nucleon off mass shell and the residual system on mass shell, are effectively 
in the same category as they have to compensate the violation of the momentum sum rule by 
introducing an additional (presumably mesonic) component in the nucleus wave function. 
These models typically predict an enhancement of the antiquark u + d distribution in nuclei 
of the order of 10-20% at x ~ 0.1 , while according to the data [3j qa/qn < 1 for this 
x-range (the absolute accuracy of these measurements is of the order of 1% ). Accounting 
for the standard model effects and accordingly reducing r/^r by a factor ~ 2 would still leave 
a gross contradiction with the Drell-Yan data. 

Another class of nuclear models is mean field models where a nucleon is treated as 
moving in the self consistent field of other nucleons with a deformation of a bound nucleon 
(nucleon swelling, ...) independent from nucleon momentum (for a review and references 
see |37j). These models ignore the experimental observation of SRCs in nuclei, and they 
do not provide an explanation for the similarity of the A-dependence of the hadronic EMC 
effect to the A-dependence of the two nucleon SRCs. The similarity of the Q^-dependence 
of the nucleon magnetic form factors of the low momentum bound nucleons to that of the 
free nucleons is not reconciled in these models with the large hadronic EMC effect. The 
fast onset of the hadronic EMC effect with x for x > 0.5 and lack of the significant effect 
at smaller x were not predicted in these models. 

In a number of models it was assumed that nonnucleonic configurations - six quark 
configurations, A isobars, etc - are present in nuclei with significant probability. If the 
SRCs were an incoherent superposition of nucleonic and non-nucleonic configurations one 
would not be able to generate an EMC effect larger than 20% of the probability of the 
SRCs, which is < 5%. This is clearly insufficient to explain the strength of the EMC effect 
of the order of 15% for x > 0.6 and A > 40. If one would assume that there exist both 
nucleonic SRC and exotic configurations with a probability > 20%, necessary to fit the 
EMC effect, the number of nucleons below the Fermi surface would drop below 60% which 
is hardly consistent with the current experience of the nuclear physics. 

To summarize, the common perception that there exist plenty of successful models of the 
EMC effect mostly arises from treating the EMC effect as an isolated phenomenon ignoring 
the totality of the constraints which were obtained from the studies of the hard nuclear 
phenomena in the past 20 years. At the moment there seems to be no viable alternative 
to scenarios where the EMC effect is associated with modifications of rare quark-gluon 
configurations selected by hard probe, and these modifications increasing with the bound 
nucleon momentum. 
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9 How nuclear medium modifies nucleon breathing 

9.1 Deformations of bound nucleon wave function within QCD 

A solution proposed in [33] is the mechanism where hard scattering off quarks selects at 
X > 0.5 rare quark- gluon configurations in bound nucleons and that the deformation of the 
bound nucleon WF is enhanced for such rare configurations. The mechanism is based on 
two fundamental and well established properties of QCD. 

Any characteristics of a composite system should fluctuate and depend on the process. 
The text book example in QED (the abelian gauge theory) is the hydrogen atom where 
the calculation of moments of the hydrogen radius finds that (< r") = J d^rr^ip'^{r) differs 
from (r)". An example of observed fluctuations in QCD is the significant difference be- 
tween the electric and axial radii of the proton: (^g ^V =0.85 fm, {r^xiai)^ ~ 0-65 fm. 
The fluctuations of strength of interaction also play a key role in the explanation of the 
phenomenon of high energy inelastic diffraction. One can prove that in QCD as in QED 
interaction of a hadron in a small size configuration with a hadron target is much weaker 
than in an average configuration. 

Important difference between nonabelian (QCD) and abelian (QED) gauge theories is 
that in nonabelian case fluctuations of wave function of a hadron include fluctuations of 
quark-gluon content of a nucleon. The strengths of the instant interaction between hadrons 
differs for different quark-gluon configurations. This is important for hard processes which 
we discuss in the next subsection. In the low energy processes instant interaction is averaged 
out and cannot reveal itself. On the contrary it follows from the QCD factorization theorem 
that a hard probe selects particular instantaneous quark-gluon configuration in a hadron 
target. Another feature of QCD is that the number of constituents is decreasing with 
decrease of the overall size. (Last property reveals itself in the x dependences of parton 
distributions, in the special high energy processes, etc) 

Difference in the fluctuations of the wave function of a bound and free nucleon leads to 
the deformation of its wave function. This phenomenon is well known in atomic physics 
from the application of the variational principle to the hydrogen molecule, H2. At large 
interproton distances the main effect is swelling of the hydrogen WFs while at small inter- 
proton distances the overall size of the hydrogen atoms is reduced [38j . 

The use of the variational principle allows us to understand the qualitative trend in 
QCD as well - probabilities of the quark-gluon configurations within a bound nucleon 
for which attraction is weaker than in average, should be reduced, while probabilities of 
configurations for which attraction is stronger than in average should be enhanced. (This 
is a particular case of the application of the Le Chatelier's principle.) 

In [15] the expressions were derived for the reduction of the probability of the configu- 
rations in bound nucleons which interact with strength much smaller than average: 

W = l + ^, (37) 
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where {U) is the expectation value of the potential acting on a nucleon in the nucleus, 
IS.E = M]\j* — M]\f ~ 400 -^ 600Mey is characteristic energy of excitation of the nucleon. 
Hence for a heavy nucleus a maximal suppression is of the order 20%. The calculation 
of the suppression effect for the case of quark-gluon configurations for which strength of 
interaction is smaller but still comparable to {U) requires use of a specific model, see |39j . 
One can use equations of motion to derive the dependence of the suppression on the 
momentum of the bound nucleon |33j. In the lowest order in k'^ /m]\fAE neglecting term 
oc ca/^E one finds 

6{k) = 1 - 2k^/mNAE. (38) 



Eq ]38| in cluding the binding term can be written in a compact form if one substitutes 
k'^ in Eqj38|by 

Am^ = m% - {pA - Precf- (39) 

where prec is the four momentum of the A-1 nucleon recoil system in the process where a 
nucleon was removed from the nucleus. The possibility to rewrite Eqj37| in this form was 
first pointed out in |40| . A detailed analysis which demonstrated validity of the formula 
for the case of generic final state of the nucleus was performed in |25j . Note here that 



Eq, 38 naturally satisfies the requirement that the scattering amplitude when continued to 
the pole where Am^ = 0, should coincide with the on-shell amplitude. Hence the linear 
dependence on off-shellness for small Am^ should be valid for wide range of bound nucleon 
deformation effects. In particular such a pattern is consistent with the measurements of 
the deviation of the ratio Q^^^^-d- / Qbound fj-Qj^ ^j^g fj.gg value which was studied in the JLab 
experiment 



9.2 Implications for nuclear pdfs at large x 

The studies of the x-dependence of nucleon pdfs find a very different dependences on large x: 
xVn{x, Q"^) cx (1 — x)^, xq]\f{x, Q"^) oc [l — xY and xGn{x, Q"^) oc (1 — x)^. This dependence 
indicates that hard interaction with nucleon at x > 0.6 selects quark - gluon configurations 
in which antiquarks in the nucleon and in particular the meson field is suppressed. The 
suppression of the interaction of nucleons in such configurations can be demonstrated in 
the perturbative QCD and on the level of emission of mesons as suppression of the pion 
emission by a small size configurations in a nucleon (see analysis in Appendix D of |15j). 

At the same time the analyses of low energy nucleon - nucleon interaction suggest that 
the dominant contribution to the SRC for k < 600 MeV/c originates from the I = 0,S = 1 
channel where pion exchange plays a very important role (see e.g. |42j ) . 

Combining this observation with the arguments of the previous subsection we conclude 
that probability of configurations with large x partons should be suppressed in the bound 
nucleons with the amount of the suppression comparable to the one given by Eq|37| which 
is of correct magnitude to explain the EMC effect for x ~ 0.6 -r- 0.7. 
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In the limit when one can average over the Fermi motion of nucleons (x < 0.6) we find 
for the A-dependence of {6): 

l-{6) = {2k'^/AE)o^TA/mNOca2{A) (40) 

Hence we conclude that for the discussed kinematics the hadronic EMC effect should be 
proportional to 02 (^) with a good accuracy which is indeed the case, see the discussion in 
section 6. 

9.3 Zooming on the nucleon polarization 

The critical test of the discussed picture would be a study of the dependence of the nucleon 
deformation on the nucleon momentum which for a large range of nucleon momenta should 
be proportional to the square of the momentum of the struck nucleon. 

The key experiment to test the discussed interpretation would be a measurement of 
the tagged EMC effect [33], the process e+'^H^e + N + X with detection of the recoil 
nucleon with light cone fraction a > 1. In this process one can measure the bound nucleon 
structure function F|^"^(x/a). We expect that 

1 - F^r''{x/a,Q^)/F2N{x/a,Q^) = f{x/a,Q')6m\ (41) 

with /(x/a,(5^) close to zero for x/a < 0.5 and rapidly growing at larger x/a. A fast 
dependence of the effect on x/a would make it easier to separate it from possible final 
state interaction effects. 

A priori the deformation of a bound nucleon can also depend on the angle (p between 
the momentum of the struck nucleon and the reaction axis as 



da 

_dn_ 

I da_\ 

\dn' 



l+cik,q). (42) 



Here {a) is cross section averaged over (j) and dil is the phase volume and the factor c 
characterizes non-spherical deformation. Correlation between the photon polarization and 
k direction is also possible. 

Such non-spherical polarization is well known in atomic physics. Contrary to QED de- 
tailed calculations of this effect are not possible in QCD. However, a qualitatively similar 
deformation of the bound nucleons should arise in QCD. One may expect that the defor- 
mation of bound nucleon should be maximal in the direction of radius vector between two 
nucleons of SRC. [^The (/>-dependent shape deformations of bound nucleons are averaged 
out in the EMC effect for x < 0.7. 

As we mentioned above a non-spherical nucleon deformation may be manifested in the 
polarization effect in the process e+^ He — )• e+p+^ H where the the measured asymmetry 



''We are indebted to H. Bethe(1995) and V.Gribov(1993) who draw attention of one of us (MS) to this 
effect. 
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can be interpreted in terms of the deviation of ratio Ge/Gm for a bound nucleon from the 
free value. It would be interesting to check whether the same pattern is observed for the 
scattering off the deuteron, whether the linear dependence on /^m? extends up to larger 
nucleon momenta and look for the dependence of the effect on the angle (p. Since the 
deformation is expected to be strongest along the radius vector between the nucleons, one 
may expect the deviations from the free nucleon case to be maximal for (/> '^ 0, vr. 

10 Conclusions 

Application of the baryon charge and momentum QCD sum rules allows us to prove that 
the EMC effect is a genuine nuclear effect. To extract hadronic EMC effect from data one 
should take into account the model independent effects: Fix the kinematic uncertainties 
in definition of Bjorken x for the scattering off nuclei as dictated by the parton model and 
QCD, and include the equivalent photon component of the light-cone WF of the heavy 
nucleus. Together these two effects account for about half of the EMC effect for x < 0.5. 
Consequently the data indicate that the hadronic EMC effect is a predominantly a high 
X phenomenon. The A-dependence of the hadronic EMC effect is stronger than for the 
ratio 2aAixp)/AaDixp) — 1 presented in the experimental studies and consistent with the 
A-dependence of the two nucleon correlations. The implementation of above mentioned 
effects impacts also on the nucleon pdfs extracted from the data and in particular strongly 
reduce the nuclear correction for extraction of F2n{x, Q'^) from the deuteron data. 

We also emphasized that the presence of predominantly nucleonic SRC in nuclei to- 
gether with the similar A-dependence of SRCs and the hadronic EMC effect, and lack 
of the nuclear enhancement of antiquark distribution in nuclei, put strong constraints on 
the mechanism of the EMC effect. These data pose a serious problem for the mean field 
models where SRC are absent. At the same time they are consistent with the scenario of 
the nucleon deformation occurring predominantly due to suppression of small size, large 
X configurations in the bound nucleons. It is based on the well established properties of 
bound states in QCD : presence of fluctuations with strength of the interaction and of 
quark-gluon content, and selection of the certain instantaneous quark-gluon configurations 
by hard processes. 

It is important to perform further experimental measurements of the EMC ratio for 
X > 0.6 at large Q^ which will become feasible at JLab 12. Especially interesting would be 
to study the ratio of the EMC effect in ^^Ca and ^"^Ca since in this case the standard model 
leads to the same nuclear effect. Moreover such research would provide a unique window 
on the structure of neutron rich high density nuclear matter relevant for description of the 
cores of the neutron stars. 

The critical tests of the current ideas will be possible in the processes where momentum 
of the struck nucleon is tagged. It is expected that effect should be proportional to the 
nucleon off-shellness and may also depend on the direction between the momentum of the 
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struck nucleoli and the virtual photon momentum (polarization). 

It would also be interesting to analyze the role of the violation of the isotopic invariance 
due to presence of the Coulomb field in nuclei and in particular the difference of the 
momentum fraction carried by photons in protons and neutrons {'r]'y{p) — "H-yin) ~ 0.2%) 
and its role in the extraction the Weinberg angle from neutrino data|43|. 
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